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Abstract 

Ab initio electronic structure studies of prototypical polar interfaces of 
wurtzite III-V nitrides show that large uniform electric fields exist in epi- 
taxial nitride over layers, due to the discontinuity across the interface of the 
macroscopic polarization of the constituent materials. Polarization fields for- 
bid a standard evaluation of band offsets and formation energies: using new 
techniques, we find a large forward-backward asymmetry of the offset (0.2 
eV for AIN/GaN (0001), 0.85 eV for GaN/AIN (0001)), and tiny interface 

formation energies. 
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Due to their low-symmetry crystal structure, wurtzite III-V nitrides exhibit a non-zero 
macroscopic polarization even in equilibrium (spontaneous polarization) .0 Because of the 
appreciable lattice mismatch between nitrides, and of the fact that nitride heterostructures 
are usually grown along the polar (0001) axis, the macroscopic polarization in an epitaxially 
grown nitride layer will include a piezoelectric term. Piezoelectric constantsB much larger 
than in most other semiconductors imply that small strains can produce unusually large 
polarizations in Ill-nitrides. Also, spontaneous and piezoelectric polarizations are compa- 
rable in magnitude.! Therefore, a major influence of polarization on interface and device 
properties should be anticipated. 

In this paper we present a detailed first-principle density functional theory study (with 
full account of strain and polarization effects) of a prototypical strained, polar, wurtzite 
nitride interface: GaN/AIN (0001). The central results discussed below are (i) the change 
in macroscopic polarization across the heterointerfaces generates large uniform electric fields 
in the layers composing the nanostructure, and (ii) a large forward-backward band offset 
asymmetry exists, due to the effects of epitaxial strain on the bulk band structure. While 
analogous (though much smaller) fields have been previously predicted in strained super- 
lattices of zincblende compounds and in ordered III-V alloys,! III-V nitrides stand alone 
because of their unusually strong polarization,0 both spontaneous and piezoelectric. The 
presence of large polarization fields has a host of interesting consequences on device design 
which will be discussed in detail elsewhere.! 

Most investigations so far have focused on the interface band offset and its possible asym- 
metry (the offset for A1N on GaN may differ from that of GaN on A1N). It is clear that 
measurements and theoretical predictions of this basic ingredient of heterostructure design 
may be significantly influenced by macroscopic polarization and by strain effects (both di- 
rect on the band bulk structures, and indirectly through piezoelectric effects). Surprisingly, 
apart from notable exceptions,!'! the recent experimental^! and theoretical!^ literature in 
this field has either disregarded or analyzed in insufficient detail the effects of macroscopic 
bulk polarization on interface electronic structure. In particular, theoretical work so far 



mostly dealt with zincblendcffl or artificially lattice-matched wurtziteEl interfaces, thus im- 
plicitly disregarding part of the key ingredients of the problem (polarization, or strain, or a 
combination of both). 

Technical details of the local-density-functional ultrasoft-pseudopotential0 plane-wave 
technique and of the theory of polarization^ employed here are reported in recent 
papers.BfflEl Results on bulk lattice parameters,0 dielectric^! and piezoelectricB constants, 
and spontaneous polarization have also been reported previously Technicalities specific to 
interface calculations will be reported elsewhere.il Here we only mention that we accurately 
reproduced previously reported studiesi'0 for GaN/AIN (111) interfaces, and that our re- 
sults for GaN/AIN (0001) are in good agreement (where they can be compared) with similar 
calculations by a different group.ll 

Here we study (GaN) n /(AlN) m (0001) superlattices such that internal fields do not cause 
metallization and at the same time the repeated interfaces are fully decoupled (n=m=4). 
Polarization effects on arbitrary nitride quantum structures will be discussed in Ref. [| We 
impose to the superlattice the in-plane lattice constant of either GaN or A1N in order to 
simulate the epitaxial relation of an hetero-overlayer on either a GaN or an A1N substrate. 
The axial lattice parameter and internal parameters of the epitaxial material are optimized 
at the imposed subtrate in-plane lattice parameter. 

We evaluate the valence band offset by splitting it conventionally^ into the difference 
AE V of the bulk valence band energies for the two bulks, and the interface potential line-up 
AV. The latter is generally just a jump in potential across the interface from one constant 
value to another. Our first result is that this is not the case at polar nitride interfaces, and 
that the line-up cannot be obtained in a conventional fashion. Indeed, consider Fig. [I], which 
show the macroscopic averag of the total charge density, and the ensuing electrostatic 
potential, of a GaN-matched GaN/AIN (0001) superlattice. The foremost unusual feature 
is of course the presence, in the bulk-like regions between the interfaces, of very large (~ 
10 9 V/m) uniform electric fields generated by the different charge distributions at the two 
interfaces (the density vanishes far from the interfaces, which indicates that the bulk-like 



regime is reached in our simulation). 

The main consequences are: (a) the difference between the bulk values of the electro- 
static potential at the two sides of the interface cannot be defined unambiguously,!^ as it 
will depend on the choice of the interface position or of the center of the bulk-like region, 
which are of course ill-defined; (b) because of energy contributions due to strain and the 
electrostatic field, the formation energy cannot be extracted as a straightforward difference 
between total energies and chemical potentials. In addition (c) the origin of the interface 
charge asymmetry must be identified; we will show that an interface charge accumulation 
takes place because of the discontinuity of the macroscopic (spontaneous and piezoelectric) 
polarization across the interface. 

We now show that the determination of the potential line-up (point (a)) and the identifi- 
cation of the sources of the uniform fields (i.e., charge asymmetry, point (c)) can be obtained 
via a multipole decomposition of the macroscopically-averaged interface charge density. The 
latter contains multipoles of all order, which in one- dimensional space are its moments. We 
are interested in the constant potential drop across the interface: this is uniquely determined 
by the interface dipole^ We are also interested in understanding the V-shaped superlattice 
potential: these are of course generated by the interface monopole^ All higher multipoles do 
not generate any potential jumps or uniform fields, but only minor potential bumps at the 
interface, symmetric and antisymmetric for even and odd multipoles respectively. Therefore, 
in practice, to extract the effects of monopoles and dipoles, we simply need to decompose 
the total macroscopically-averaged charge density n into two components comprising, re- 
spectively, all its even and odd multipoles. 

For the sake of clarity, we name the odd and even components respectively the dipole 
density Edi p , and the monopole density n mono . This is admissible since these densities produce 
all the effects of dipolar and monopolar charges relevant to our problem, plus other minor 
effects related to higher multipoles (irrelevant for our purposes). 

Unfortunately, such a decomposition can be done in an infinite number of ways. Our pro- 
cedure to obtain n mono is to fold the density with respect to a mirror plane placed at a point 
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zq roughly halfway between two adjacent interfaces, and then perform an antisymmetric 
combination of the two charge distributions thus superimposed, i.e. 

™mono(-Z - Zq) = - {n(z - Zq) ~ Tl(z - z)) , (1) 

where z is the position of the folding plane. The dipole density ndi P is defined as the 
difference between the full density and the monopole term, or (which amounts to the same) 
the symmetric combination of the two superimposed densities 

n d i P (z - Zq) = ^ [n(z - Zq) + n(z - z)} . (2) 

The key point of this procedure is of course the choice of the position zq for the folding 
plane, which implicitly selects one specific realization of the decomposition. Our criterion 
for choosing z is that the norm of the monopole component, 

S(z ) = J \n moao (z - z )\ 2 dz , (3) 

should be minimized. This choice produces (a) a dipole distribution which deviates mini- 
mally in a least-squares sense from the total density; (b) a n mon o optimally localized at the 
interface; (c) a position for the folding plane that coincides with the intuitively appealing 
idea of midpoint between adjacent interfaces. 

The monopole and dipole distributions obtained by the above decomposition are shown 
in Fig. H] for a typical case. The dipole is related to a jump in potential across the interface, 
and it allows the direct determination of the line-up potential, and therefore of the band 
offset. In turn, the interface monopole can be further analyzed to ascertain its physical 
origin. 

Let us first present the valence-band offset of the GaN/AIN (0001) interface. The offset 
is of type I. As reported in Table |, we obtain A_E„=0.20 eV for A1N lattice matched 
to GaN, and A£'„=0.85 eV for GaN matched to A1N. We thus confirm the existence of 
a large forward-backward asymmetry (0.65 eV) of the offset. This asymmetry originates 
mostly from strain-induced band edge shifts in the bulk band structures (contributing 0.47 

5 



eV), with moderate contributions from the lineup term (0.18 eV). The large band offset 
asymmetry is thus mostly to be ascribed to the different nature of the valence band edge 
in GaN and A1N. Indeed, the A1N (GaN) band edge is a singlet (doublet) formed by the 
hybridization along the c-axis (in the a-plane) of N 2s orbitals with Al p z (Ga p xy ) states, so 
that biaxial compression pushes the edges upward in GaN and downward in A1N. We note 
that our findings are semiquantitatively in agreement with those reported by Nardelli et a/.i 
for zincblende (001) interfaces, namely 0.44 eV for A1N on GaN, and 0.73 eV for GaN on 
A1N. 

Let us now turn to the interface monopole. The dipole is understood to belli a response 
to the electrostatic perturbation induced by interface formation (for the present system, 
in which this effect is adulterated by monopole contributions, our decomposition gives the 
best approximation to this response). On the other hand, the monopole may be expected 
to be the difference in macroscopic polarization between the constituents of the junction: 
indeed, according to Poisson's equation, a polarization discontinuity at the interface be- 
tween two different media produces an interface charge accumulation. In particular, in a 
superlattice made of alternating layers of materials A and B of respective thicknesses I a and 
Ib and dielectric constants Ea and eb, the areal charge density at the interface is directly 
connected@0il with the transverse bulk polarizations Pj and Pg of the interfaced materials 
by 

<Tint = (Pj - Pi) Ha + Ib)/(Ia£b + l B e A ), (4) 

where we have assumed conventionally that Pg (Pj) is the transverse polarization on the 
right (left) side of the interfaced This relation allows an independent prediction of what the 
polarization-induced interface monopole should be, which can be compared with the minimal 
monopole calculated for the actual interface. Fortunately, the transverse polarization P T 
of the nitrides can be computed accurately^ by means of the geometric quantum phase 
approach@ in an arbitrary strain state, for instance for the epitaxially strained overlayer 
material. The dielectric constants (static or electronic) can also be evaluated independently 



using a recently developed technique.El 

For the unrelaxed structure (clamped ions), the electronic dielectric constant should be 
used in Eq. [|, as appropriate to purely electronic screening. In the real system, however, 
the electric field induces a lattice distortion which extends over the whole slab, i.e. a long- 
wavelength optical phonon gets frozen-in: it is then appropriate to use in Eq. |] the static 
dielectric constant as calculated in our previous work.0 In Table | we report the actual 
interface charge density obtained via the multipole decomposition, and the value cr i ^ P ' ) 
obtained from Eq. f|, for both the ideal and the relaxed superlattice. The excellent agreement 
of the pairs of independently determined values confirms indeed the identification of the 
interface charge with a polarization charge. 

A final important issue is the evaluation of the interface formation energy. For the 
present system it is impossible to build a superlattice with equivalent interfaces, so that 
a total energy calculation can only provide an average interface formation energy; this is 
hardly a severe problem, as the two interfaces are very similar.0 An additional problem is 
that the superlattice total energy contains elastic and electrostatic energy contributions due 
to, respectively, lattice mismatch and polarization fields. Clearly, these contributions are 
extensive, i.e. they depend on the overlayer thickness when referring formation energies to 
unit area. In analogy to surface energies, we write the total energy per superlattice unit cell 
as 

£ t s > x ) = 2Ef + £ n x (^ + e + V X ) (5) 
x 

where /i x are the total bulk energies per Ga-N or Al-N pair (in the appropriately strained 
geometries), £ x are the elastic energies and r] x the electrostatic energies stored in the (pos- 
sibly) strained bulks under the polarization field, and n x is the number of atom pairs of 
type X (GaN or A1N). 

In the present case of a strained low-symmetry system, an exact numerical equivalence of 
bulk and interface (in particular between k-point meshes) cannot be achieved, and the use of 
fi, £ and rj evaluated from separate bulk calculations might lead to inaccuracies. A solution 



to this issue, as in the case of surfaces jHil is to recognize that Ef^ depends linearly on n x , 
so that E* can be extracted as the intercept of the linear E tot vs. nx relation, i.e. from a 
series of total energy calculations for superlattices of different lenghts (whereby equivalent 
k-point sets are easily obtained). 

In Table [TJ we list the formation energies for the ideal and relaxed interfaces obtained 
by linear extrapolation. The same Table reports bulk values of the elastic and electrostatic 
energy, the former obtained as total energy difference with the unstrained lattice, and the 
latter as 

V X = \e x n*E* (6) 

with Q x is the bulk cell volume, E the modulus of the electrostatic field, and ex the static 
dielectric constant of material X (the static dielectric constant implicitly accounts for the 
field- lattice coupling!!!) . The strain energy is much larger than the interface energy and the 
electrostatic energy, even for modest thicknesses. Assuming an order of magnitude for the 
dislocation core formation energy of ~ 0.5 eV,il we see that the formation of such strain- 
related defects should start at typical thicknesses of ~ 20 A. A comparable electrostatic 
energy would be stored in (perfect) layers of thickness in the order of 500 A, and is therefore 
irrelevant to the layer's stability, since metallization or screening effects set in at much 
smaller thicknessesJl Thus, it can be safely stated that (0001) nitride interfaces are abrupt, 
and that the electrostatic energy should not prevent their stability, as it may in heterovalent 
systems such as ZnSe/GaAs. 

In summary, our study of AlN/GaN(0001) interfaces has revealed the presence of large 
uniform electrostatic fields which we demonstrated to originate from the macroscopic polar- 
ization of the junction constituents. We have also indicated ways of extracting band offsets 
and formation energies, for which conventional definitions are useless in the present situa- 
tion. We found a sizable forward-backward band offset asymmetry, tiny interface formation 
energies, and large epitaxial strain energies. 
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TABLES 



Substrate -> GaN A1N 

(0.29) 0.85 (1.00) 

(0.029) 0.011 (0.022) 

(0.028) 0.011 (0.022) 

TABLE I. Valence-band offset AE V (eV) and mono-pole charge cr int (C/m 2 ) at AIN/GaN (0001) 
for different epitaxial matching conditions, and fully relaxed super lattices (in parenthesis: unre- 
laxed case). 



Substrate [ 


Ef 


V AIN 


^GaN 




^GaN 


GaN 


3.9 


5.6 


9.7 


179 




A1N 


0.4 


10.9 


6.3 




155 



TABLE II. Formation, electrostatic, and elastic energy for an AIN/GaN superlattice for differ- 
ent substrate choices (meV/cell or unit area). 



(SL) 
int 



(AP) 



int 



0.20 

0.014 

0.014 
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FIGURES 

FIG. 1. Total (electronic plus ionic) density and ensuing electrostatic potential (in Hartree) for 
an AIN/GaN super lattice matched-lattice to GaN. The magnitude of the fields in the bulk regions 
is ~ 10 9 V/m. 

FIG. 2. Full density (dash-dotted), and monopole (solid) and dipole (dashed) components for 
the superlattice of Fig.|l]. 
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